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Abstract 

Understanding glasses is considered to be one of the most fundamental problems in statistical 
physics. A theoretical approach to unravel their universal properties is to consider the validity of 
equilibrium concepts such as temperature and thermalization in these out-of-equilibrium systems. 
Here we investigate the autocorrelation and response function to monitor the aging of a colloidal 
glass. At equilibrium, all the observables are stationary while in the out-of-equilibrium glassy state 
they have an explicit dependence on the age of the system. We find that the transport coefficients 
scale with the aging-time as a power-law, a signature of the slow relaxation. Nevertheless, our 
analysis reveals that the glassy system has thermalized at a constant temperature independent 
of the age and larger than the bath, reflecting the structural rearrangements of cage-dynamics. 
Furthermore, a universal scaling law is found to describe the global and local fluctuations of the 
observables. 
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Introduction. — Increasing the volume fraction of a colloidal system slows down the 
Brownian dynamics of its constitutive particles, implyinga limiting density, </)g, above which 
the system can no longer be equilibrated with its bath Hence, the thermal system falls 
out of equilibrium on the time scale of the experiment and thus undergoes a glass transition 
j^. Even above (pg the particles continue to relax, but the nature of the relaxation is very 
different to that in equilibrium. This phenomenon of a structural slow evolution beyond the 
glassy state is known as "aging" P|. The system is no longer stationary and the relaxation 
time is found to increase with the age of the system, t^, as measured from the time of sample 
preparation. 

This picture applies not only to structural glasses such as colloids, silica and polymer 



astic manifolds in quenched disor- 
Theories originally developed 



melts, but also to spin-glasses, ferromagnetic coarsening, e 
der and jammed matter such as grains and emulsions p 
in the field of spin-glasses attempt to develop a common framework for the understanding 
of aging. For example, the structural glass and spin-glass transitions have been coupled by 
the low temperature extension of the Mode-Coupling Theory (MCT) jj, 0, Q]- More gen- 
erally, this approach is related to analogous ideas developed in the field of granular matter 



such as compactivity 
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to the energy landscape of glasses 
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and the inherent structure formalisms 



0, 



adapted 



One of the important features of this scenario is a separation of time-scales where the 
observables are equilibrated at different temperatures, even though the system is far from 
equilibrium. While theoretical results have flourished, the difficulties in the experimental 
testing of the fundamental predictions of the theories have hampered the development of 
an understanding of aging in glasses Q|- Experiments so far have shown 

conflictin g re sults which are usually masked by large intermittent fluctuations in the ob- 



a behavior which seems beyond the current theoretical formalisms 
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servables 

the other hand, some numerical results are more favorable 

0,0, Q- 

Furthermore, the 

concept of temperature has been shown to be useful to describe other far form equilibrium 
systems such as non-thermal granular materials 0, ^ , U , 3] • 

Here we use a model glass which is one of the simplest systems undergoing a glass transi- 
tion: a colloidal glass of micrometer size particles, where the interactions between particles 
can be approximated as hard core potentials 1, 23, 2^. The system is index matched to al- 
low the visualization of tracer particles in the microscope 2^. Owing to the simplicity of the 
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system, we are able to follow the trajectories of magnetic tracers embedded in the colloidal 
sample and use this information as an ideal "thermometer" to measure the temperature for 
the different modes of relaxation. In turn, we measure the autocorrelation function of the 
displacements and the integrated response to an external magnetic field as an indicator of 
the dynamics via a Fluctuation-Dissipation Theorem (FDT). For this system we show that, 
even though the diffusivities and mobilities of the tracers scale with the age of the system, 
there exists an effective temperature which is uniquely defined and remains constant inde- 
pendent of the age. This effective temperature is larger than the bath temperature and 
controls the slow relaxation of the system, as if the system were at "equilibrium" . We find a 
scaling behavior with the waiting time, which describes in a unified way not only the global 
but also the local fluctuations of the correlations and responses as well as the cage dynamics 
in the system. 

Correlations and responses. — Our experiments use a colloidal suspension consisting 
of a mixture of poly-(methylmethacrylate) (PMMA) sterically stabilized colloidal particles 
(radius = 1.5 /im, density pp = 1.19 g/crn?, polydispersity ~ 14%) plus a small frac- 
tion of superparamagnetic beads (radius = 1-6 /im and density pm = 1-3 g/cm?, from 
Dynal Biotech Inc.) as the tracers (More details in the Method section). In order to 
investigate the dynamical properties of the aging regime, we first consider the autocorre- 
lation function as the mean square displacement (MSD) averaged over 82 tracer particles, 
C(t,tw) = {Ax'^(t,tw)) /2 = {[x{tw + At,tw) — x{tw,tw)Y) /2, at a given observation time, 
t = tw + At, after the sample has been aging for as measured from the end of the 
stirring process (see Appendix [HI). Then, we measure the integrated response function 
(by adding the external magnetic force, F) given by the average position of the tracers, 
X(t,tui) = {x{tw + At,tw) - x{tw,tui))/F. 

Analytical extensions of the MOT for supercooled liquids to the low temperature regime of 

E lasses allow for the interpretation of the aging of global correlation and response functions 
], y. In these frameworks, the evolution of C(t,tw) and xit^tw) are separated into a 
stationary part (short time) and an aging part (long time): C{t, t^) = Cst(t — t^) + Cag(t, t^) 
and xit^tw) = Xstit — tiu) + Xagit, t^), where we have included the explicit dependence on 
in the aging part. 

This result can be rationalized in terms of the so-called "cage dynamics" . As the density 
of the system increases the particles are trapped in cages. The motion inside the cage is still 
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equilibrated at the bath temperature and is determined by the Gibbs distribution of states. 
This dynamics gives rise to the stationary part of the response and correlation functions, 
which satisfy the usual equilibrium relations such as the FDT. However, the correlation does 
not decay to zero but remains constant since particles are trapped in cages for a long time. 
Thermal activated motions lead to a second structural relaxation which is responsible for the 
aging part of the dynamics. In this regime the system is off equilibrium and correlations and 
responses depend not only on the time of observation t but also on the waiting time, tw Slow, 
non-exponential relaxation ensues and a strong violation of the FDT is expected. However, 
the interesting result is that this breakdown leads to a new definition of temperature for the 
slow modes which has been proposed to be the starting point of a unifying description of 
aging in glassy systems 

Figure^ shows {Ax'^{t,tw)) as a function of At at a fixed = 100s calculated for the 
three colloidal samples at (pc < 'Pg < 'Pa < 'Pb, and Fig. [Hd shows the age dependence for 
the glassy sample A. Sample B shows t^j dependence similar to sample A, while sample C, 
being at equilibrium, is independent of t^, i.e. it is stationary (see Appendix lE| . 

The cage dynamics is evidenced by the plateau observed in the MSD in the two glassy 

2511 to be observed 



samples. The rattling of particles inside the cages is too fast (~ 10~ 
with our visualization capabilities. Since we focus mainly on the long relaxation time regime, 
the stationary parts, Cst and Xst, are negligible compared with Cag and Xag, in the following 
we concentrate only on the aging part of the observables and drop the subscript ag: C = Cag 
and X = Xag- The tracers' motion is confined by the cage, which persists for a time of the 
order of the relaxation time r(t^). As expected for an aging system, this relaxation time 
increases with as observed in Fig. ^p. For longer times. At > r(t^), structural rearrange- 
ments lead to a second increase of the MSD, defining a diffusion regime characterized by a 
diffusion constant, D{t^,), which depends on the waiting time. We find an asymptotic form: 

{[x{t^ + At,t^)-x{t^,t^)f)r^2D{t^)At, foi At >r{t^). (1) 

Figure shows the average displacement of the magnetic beads under the external force 
as a function of time At, for various aging times, t^, in sample A. The magnetic force is set 
as small as possible to observe the linear response regime, F = 1.7 x 10~^^A^ (see Appendix 
IF|) . Contrary to the behavior of the MSD, the integrated response function does not display 
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the plateau characteristic of the cage effect. The data can be fitted to: 

{x{t^ + At, t^) - x{t^, t^)) ~ M{t^)FAt, (2) 

where M(t^) is the mobihty of the tracers which is again waiting time dependent as seen in 
the figure. 

To investigate the nature of the scahng behavior of the aging regime, we study the 
dynamical behavior with respect to t^. Figure |21 shows both the diffusivity and mobility as 
a function of for sample A. Both D{t^) and M{t^) decrease with signaling the slowing 
down in the dynamics. More importantly, they decrease according to a power-law with the 
same exponent for both quantities: 

D(t^) oc tw''^ and M(t^) oc t^"''', (3) 

with 7 = 0.32 ±0.08. This result is consistent with previous work in a similar agin g co lloidal 
system, where it was found that the MSD (Ax^) has a power law decay with j26[. The 
fact that all the quantities scale as power laws indicates that the aging regime lasts for a 
very long time, perhaps without ever equilibrating. 

Effective temperature. — The same power law decay of the diffusivity and mobility 
implies that the system has thermalized at a constant effective temperature Teg independent 
of tuj. This temperature is given by an extension of the Stokes-Einstein relation or FDT to 
out-of-equilibrium systems. Even though both D and M depend on the age of the system, 
their ratio is constant: 

= riT^fr = ^690 ± 100)K. (4) 

The inset of Fig. |21 plots the function of t^. We obtain Tcs ~ 690K which is 

more than double the ambient temperature of 297K. For very large ty^, T^g shows large 
fluctuations that are mainly due to the larger statistical error (due to the limited number 
of tracers) in obtaining the mobilities and diffusivities in the large At and regime. It 
remains a question whether the long waiting time regime may show interrupted aging. 

Equation (0]) is easy to understand when the system is at equilibrium: we extract energy 
from many identical tracers located in distant regions of the colloidal system and transfer 
it to the thermometer system. The thermometer receives work from the diffusive motion of 
the tracers, and it dissipates energy through the viscosity of the system. These two opposing 
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effects malce the thermometer stabihze at a temperature guaranteed by the Einstein relation. 
Naturally, we have applied the diffusion-mobility calculations to dilute sample C and find 
that it is equilibrated at the bath temperature (see Appendix . On the other hand, the 
colloidal sample is aging out of equilibrium. Nevertheless, the fact that the ratio of diffusion 
to mobility yields a constant temperature can be taken as an indication that the long time 
behavior of the system has thermalized at a larger effective temperature Tcs ~ 690K. 

Although it may seem counterintuitive that the slow relaxation at long times corresponds 
to a temperature that is actually higher than the equilibrium bath temperature, there is an 
interesting physical picture that rationalizes this observation. One can think of the energy 
landscape of configurations of the colloidal glass being explored less frequently, yet the 
amplitude of the jumps between basins corresponds to "hotter" explorations of a boarder 
distribution of energy states. While this mechanism violates the usual relations between 
particle motion and temperature, it gives rise to the effective temperature measured in our 
experiments. 

Scaling ansatz for the global correlations and responses. — Further insight into 
the understanding of the slow relaxation can be obtained from the study of the universal 
dynamic scaling of the observables with tw Based on spin-glass models, different scaling 



scenarios have been proposed [2, |27| for correlation and response functions. Our analysis 



indicates that the observables can be described as 



C{t^ + At,t^) = {Ax')/2 = C/d(-^), (5a) 

(Ax) At 
XK + At, t^) = ^—^ = t-"/A,(— ), (5b) 

where fn and are two universal functions and a and (3 are the aging exponents. Evidence 
for the validity of these scaling laws is provided in Fig. |21 where the data of the correlation 
function and the integrated response function collapse onto a master curve when plotted 
as f^Cityj + At.tu]) and t^x(^«« + ^i^iw) versus At/t^^. By minimizing the cx^ value of 
the difference between the master curve and the data (see Appendix EI) we find that the 
best data collapse is obtained for the following aging exponents: a + (3 = 0.34 ± 0.05 and 
(5 = 0.48 ±0.05. We find (Fig. ^ that the scaling functions satisfy the following asymptotic 
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behavior: 



y 



,0.3 



y < 1, 



foiy) ~ { 



(6a) 



y 



y > 1, 



(6b) 



in agreement with the fact that the motion of the particles is diffusive at long times, (Aa;^) ~ 
At, and the existence of a well-defined mobility, respectively. Therefore at long times, both 
the correlation and response functions display the same power law decay: 



The result C{t^ + At,t^) ~ x{tw + At,t 

~ ^w^"^^ confirms our previous result, Eq. 
dSl), D{t^) ~ M{t^) ~ This is in further agreement with the finding that T^s 

is independent of the age of the system, t^. For short times, the MSD scaling function 
crosses over to a sub-diffusive behavior of the particles. We obtain, {x'^ity, + At,t^)) ~ 
t~'^{At/t^)^'^ = t~°'°°^At°''^. Since the trapping time corresponds to the size of the cages 
denoted by q{tyj), we can determine the cage dependence on as qit^) ~ The 
resulting exponent is so small that we can say that the cages are not evolving with the 
waiting time, within experimental uncertainty. Furthermore, the scaling ansatz of Eq. ^ 
indicates that the relaxation time of the cages scales as r(t^) ~ since this is the time 
when the sub-diffusive behavior crosses-over to the long time diffusive regime. 

From Fig. El we see that the asymptotic At-linear regime appears when the reduced 
variable iStjtJ^ is larger than 10, t^tjtj^ > 10. The time separation between At and ty^ can 
be determined using the cut-off: At/tyf-'^^ ~ 10. For smaller times, At/ty,^ < 10, we obtain 
a subdiffusion regime (with exponent 0.3), which is characteristic of the cage dynamics. For 
longer times, AtjiJ^ > 10, we observe the asymptotic At-linear regime where the diffusivity 
is calculated. The measurement of MSD in Fig. 1 extends up to = 6000s. For this t^, the 
separation of time scales appears when At > 10 * 6000°'^^ ~ lO^s. Our measurements for 
MSD extend to lO^s, ensuring a separation of time scales even for this longest waiting time. 
For the smaller ty^ considered in the calculation of Tgff , the separation of time scales is even 
more pronounced. For instance, for a typical tyj = 1000s where the T^tf is calculated, the 
separation of time scales occurs at 10 * 1000°'^^ ^ 275s, again ensuring a well defined long 
time asymptotic behavior for At = lO^s. 



C{ty,+At,ty,)r~.t-(- + '')At, 

X{ty,+At,ty,)r~.t-(-+P)At. 



(7a) 

(7b) 
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Although the existence of an effective temperature can be rationahzed using theoret- 
ical frameworks of disordered spin-glass models ^, we find that the scaling forms of 
the correlations and responses are not consistent with such models. Based on invariance 
properties under time reparametrization, spin-glass models predict a general scaling form 
Cag(t,t^) = Ci^g{h(t)/h{t^)), where h{t) is a generic monotonic function We find that 
the scaling of our observables from Eqs. (j3aj) . ()5b|) cannot be collapsed with the ratio 
h{t)/h{tw)- The scaling with h(t)lh(t,„) is expected for system in which the correlation 
function saturates at long times 2M- other hand, our system is diffusive, and the 

studied correlation function is not bounded. Indeed, similar scaling as in our system has 
been found in the ag ing dynamics of anther unbounded system: an elastic manifold in a 



disordered media [29[. The suggestion is that this problem and the particle diffusing in a 
colloidal glass may belong to the same universality class. Furthermore our results can be 
interpreted in terms of the droplet picture of the aging of spin glass, where the growth of 
the dynamical heterogeneities control the aging. 

Local fluctuations of autocorrelations and responses. — Previous work has re- 
vealed the existence of dynamical heterogeneities, associated with the cooperative motion of 
the particles, as a precursor to the glass transition as well as in the glassy state 
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2J,|2al3Q|. 



Instead of the average global quantities studied above, the existence of dynamical hetero- 
geneities requires a microscopic insight into the structure of the glassy. Earlier studies 
focused mainly on probability distributions of the particles displacement near the glass 
transition. More recent analytical work in spin glasses shows that the probability distri- 
bution function (PDF) of the local correlation P{C) and the local integrated response -P(x) 
could reveal essential features of the dynamical heterogeneities. 

Here we perform a systematic study of P{C) and -P(x) in sample A, and the resulting 
PDFs are shown in Fig. |3] The scaling ansatz of Eq. (0) implies that P{C) and P{x) should 
be collapsed by rescaling the time At by and the local fluctuations by (see Appendix 
IHfor more details). Indeed, this scaling ansatz provides the correct collapse of all the local 
fluctuations captured by the PDFs, as shown in Fig. E^, andEh for P{C) and in Fig. HJi 
for P(x). 

The PDF of the autocorrelation function displays a universal behavior following a mod- 
ified power-law t~"P(C) oc (t^C + Cq)"'^, where Cq and A only depend on the time ratio 
At/t^. For the smaller values of C (C < t~"Co), the existence of a fiat plateau in P{C) 
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indicates that the tracers are confined in the cage. For larger values of C, the salient feature 
of the PDF is the very broad character of the distribution, with an asymptotic behavior 
P{C) ~ C~^. This large deviation from a Gaussian behavior is a clear indication of the 
heterogeneous character of the dynamics. Furthermore, the exponent A decreases from 2.6 
to 1.4 with the time ratio At/t(^ ranging from 10 to 60. We notice that A = 2 corresponds 
to the crossover between the short-time and long-time regime in Fig. El where At/t£ ^ 40. 
The significance of A = 2 is seen in the integral / P{C)CdC. For A > 2 {At/tf^ < 40) the 
plateau dominates over the power law tail in the integral and the dynamics is less hetero- 
geneous. For A < 2 {At/t^ > 40) the power law tail dominates and this regime corresponds 
to the highly heterogeneous long-time regime (see Appendix QJ. 

On the contrary, P{x), shown in Fig. EJi, displays a different behavior. The fluctuations 
are more narrow and the PDF can be approximated by a Gaussian. This is consistent with 
the fact that we did not flnd cage dynamics for the global response in Fig. Moreover, 
numerical simulations of spin-glass models seem to indicate a narrower distribution as 
found here. 

We have presented experimental results on an aging colloidal glass showing a well-defined 
temperature for the slow modes of relaxation of the system. This T^s is larger than the 
bath temperature since it implies large scale structural rearrangements of the particles. In 
other words, it controls the cooperative motion of particles needed to relax the cages. The 
interesting result is that this temperature remains constant independent of the age, even 
though both the diffusivity and the mobility are age dependent. The power-law scaling 
found to describe the transport coefficients indicates the slow relaxation of the system. A 
universal scaling form is found to describe all the observables. That is, not only the global 
averages, but also the local fiuctuations. The scaling ansatz, however, cannot be described 
under present models of spin-glasses, but it is more akin to that observed in elastic manifolds 
in random environments suggesting that our system may share the same universality class. 
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Method: Experimental details 

The colloidal suspension is immersed in a solution containing 76% weight fraction of 
cyclohexylbromide and 24% cis-decalin which are chosen for their density and index of 



refraction matching capabilities 



0.57 - 0.5^ 
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2J|. For such a system the glass transition occurs at 



\. In our experiments we consider three samples at different densities 
and determine the glassy phase for the samples that display aging. The main results are 
obtained for sample A just above the glass transition (px = 0.58 ± 0.01. We also consider 
a denser sample B with 0b = 0.60 ± 0.01, although this sample is so deep in the glassy 
phase that we are not able to study the slow relaxation of the system and the dependence 
of the waiting time within the time scales of our experiments. Finally, we also consider a 
sample C below the glass transition 0c = 0.13 ± 0.01 < (pg for which we find the usual 
equilibrium relations (see below). Prior to our measurements, the samples are homogenized 
by stirring them for two hours to achieve a reproducible initial time (see Appendix El for a 
full discussion). 

We use a magnetic force as the external perturbation to generate two-dimensional motion 
of the tracers (see Appendix El) on a microscope stage following a simplified design of 
js^ . Video microscopy and computerized image analysis are used to locate the tracers in 
each image. We calculate the response and correlation functions in the x-y plane. 
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FIG.^ Autocorrelation and response functions, (a) Mean square displacement of tracers 
as a function of At for the three samples (^c* = 0.13 < (pg < (pA = 0.58 < (j)B = 0.60 at 
tyj = lOOs. The straight dash lines indicate estimates for the diffusivity D{t^). (b) Mean 
square displacement of tracers in sample A as a function of time At, for various aging times 
tu,- The dashed straight lines indicate the fitting regime to calculate the diffusivity Dit^). 
(c) Mean displacement of tracers under the magnetic force in sample A for various aging 
times tw 

FIG. 121 Diffusivity (red dash-dot line) and mobility (blue dot line) as a function of aging 
time for sample A. A straight dash line is added to guide the eyes, which shows D(ty,) ~ t^'"^ 
and M(t^) ~ with 7 = 0.32 ± 0.08. For convenience of comparison with diffusivity, 
mobility is scaled by ksTcs, with kB the Boltzmann constant and the effective temperature 
Teff = 6907^. The inset shows the effective temperature as a function of waiting time. Error 
bars are added only to some representative points for clarity. 

FIG. 13 Scaling plot of sample A for the scaled autocorrelation, tu,"C, and scaled inte- 
grated response, ksT^stw^Xi ^ function of the time ratio, At/tw^, for different waiting 
times. The black dash line is a linear fit which indicates that Tcs = 690K. The inset is a 
plot of sample B for autocorrelation, C, and integrated response, ksTcsX^ ^ function of 
At at tyj = 100s. The black dash line is a linear fit which indicates that T^s = 1600/^. 

FIG.H (a) PDF of the scaled local correlation t^C for At/tf, = 10 with At from 155s to 
650s. Solid line corresponds to the modified power-law fit t~"P(C) oc (t°C -|- Cq)"^, with 
Co = 500 and A = 2.6. (b) Same as (a), for At/t^ = 40 with At from 620s to 2600s. The 
resulting distribution can be fitted (solid line) by a modified power law with Co = 500 and 
A = 1.9. (c) Same as (a), for At/t(^ = 60 with At from 930s to 3900s. The solid line is 
the best fit, with parameters Co = 300 and A = 1.4. (d) PDF of the scaled local integrated 
response for At/t^ = 1 shows a Gaussian decay. For the convenience of the comparison with 
(a), (b), (c), we rescale the response by ksTcff, where T^s = 690K. 
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SUPPLEMENTARY MATERIALS 



APPENDIX A: DETAILS OF THE EXPERIMENTAL SET UP 

The experimental set up is shown in Fig. and Fig. Et". We use a Zeiss microscope 
with a 50x objective of numerical aperture 0.5 and a working distance of 5mm. We work 
with a field of view 194/im x 155fim by using a digital camera with the resolution of 1288 
pixels X 1032 pixels. We locate the center of each bead position with sub-pixel accuracy, by 
using image analysis. For the condensed samples A and B, we use the low frame rate 1/3 
frame/sec, while for the dilute sample C, we record the images at 1 frame/sec. The long 
working distance of the objective is necessary to allow the pole of the magnet to reach a 
position near the sample. An example of the images of the tracers obtained in sample A is 
shown in Fig. where we can see the black magnetic tracers embedded in the background 
of nearly transparent PMMA particles. An example of the trajectory in the x-y plane of a 
tracer diffusing without magnetic field is shown in Fig. Eli. We note that this particular 
tracer moves away from two cages in a time of the order of 4 hours. 

The magnetic field is produced by one coil made of 1200 turns of copper wire. We arrange 
the pole of the coil perpendicular to the vertical optical axis, and generate a field with no 
vertical component. Thus, the tracers move in the x-y plane with a slight vertical motion 
which is generated by a density mismatch between the tracers and the background PMMA 
particles. This vertical motion is very small at the high volume fraction of interest here, and 
therefore we calculate all the observables in the x-y plane. 

The magnetic force is calibrated for a given coil current by replacing the suspension with a 
mixture of 50 : 50 water-glycerol solution with a few magnetic tracers. The distance between 
the top of the magnetic pole and the vertical optical axis is always fixed, which means that 
the magnetic force at the local plane depends only on the coil current. At a given current, we 
determine the velocity of the magnetic tracers at the focal point and calculate the magnetic 
force from Stokes's law, F = GnrjamV, where rj is the viscosity of the water-glycerol solution, 
ttm is the tracer radius, and v is the observed velocity of the tracer. The uncertainty in 
the obtained force comes from: (a) the uncertainty in the coil current which is 1%, (b) the 
beads, which are not completely monodisperse in their magnetic properties, causes a 10% 
uncertainty [s^, and (c) the magnetic field is slowly decaying in the field of view, causing a 
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FIG. 5: (a) Picture of the experimental setup, (b) Schematic illustration of the setup, (c) Detail of 
an image of tracer in sample A. (d) Trajectory of a tracer in sample A showing the cage dynamics 
over 4 hours. The circle represents the size of the tracer particle, 3.2fim. 

4% uncertainty. 

APPENDIX B: SAMPLE PREPARATION AT = 

It is important to correctly determine the initial time of each measurement to reproduce 
the subsequent particle dynamics. We initialize the system by stirring the sample for two 
hours with an air bubble inside the sample (see Fig. EJd) to homogenize the whole system 
and break up any pre-existing crystalline regions. Then we place the sample on the magnetic 
stage and take images to obtain the trajectories of the tracers that appear in the field of view. 
The initial times t = ty^ = are defined at the end of the stirring. After measuring all the 
tracers appearing in the field of view, a new stirring is applied, the waiting time is reset to 
zero, and the measurements are repeated for a new set of tracers. We have analyzed the pair- 
distribution function, g{r), and two-time intensity autocorrelation function, g2{tw,tw + At), 
in order to test our rejuvenation technique. 
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2 4 6 8 10 

FIG. 6: Pair-distribution functions of (a) sample A and (b) sample B. 

First, Fig. IHl plots the pair-distribution functions of sample A and sample B right after 
the stirring procedure. We calculate the pair-distribution functions by reconstructing the 
packings from 3D confocal microscopy images of size 60/imx60/imxl5/im. We find that 
the samples do not show obvious crystallized region by directly looking either at the pair- 
distribution function or the images taken from confocal microscopy. For these measurements 
we use a Leica confocal microscope. The PMMA particles are fluorescently dyed so that 
they are ready to be observed by confocal microscopy. We load the samples sealed in a glass 
cell on the confocal microscope stage and use a Leica HCX PL APO 63x, 1.40 numerical 
aperture, oil immersion lens for 3D particle visualization in order to calculate the volume 
fraction and the pair-distribution function. 

Second, we analyze g2{tw + ^t, t^) and show that the sample has been rejuvenated by the 
stirring process. We plot g2{tw + ^t,tj^) for 3 situations: (a) after the first stirring process 
(tw = 0), (b) after the sample has aged for a long time {t^ = 4000s), (c) we then apply 
the stirring again and immediately plot the autocorrelation function for the new initial time 
{tj = 0). The plot show the lack of correlation for the two individual measurements = 
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FIG. 7: Two-time intensity autocorrelation function, g2{tw + ^tjtw), of (a) sample A and (b) 
sample B. The solid black and navy curves correspond to ty^ = Os and ty^ = 4000s, respectively. 
The dash gray and light blue curves correspond to the measurements after the second stirring with 
tw = Os and t^' = 4000s, respectively. The fact that the curves at t^, = Os and tj = Os coincide 
indicates that sample has been fully rejuvenated to its initial state. The inset is the same plot as 
the main figure for different tw ranging from Os (solid black) to 4000s (solid navy). 

and tj = 0, thus demonstrating that the sample has been rejuvenated. Technically, we 
record the temporal image sequence of sample A and sample B, and study the aging by 
calculating two images' correlation defined as: 

g2{t^n,tyj + At) = , (Bl) 

where / is the average gray-scale intensity of a small box (20 pixels x 20 pixels, PMMA 
particles' diameter is roughly equal to 20 pixels), and I{tw) and I{t) come from the images 
at different times t^j and t, respectively. We cut one image (1288 pixels x 1032 pixels) into 
many boxes, and calculate the correlation of two boxes located at the same position in the 
two images. The average (...) is taken over all the boxes in one image. Fig. [7| plots the 
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correlation function, g2{tw^ ty^ + At), of sample A and sample B, which shows aging behavior 
(see the inset). Fig. [7| shows how g2{'tw,'tw + At) calculated after two stirring processes at 
= and tj = coincide, indicating that the sample can be fully rejuvenated to its initial 
state by our stirring technique. 

APPENDIX C: DATA ANALYSIS 

The average (■ ■ ■) used to calculate the MSD (Ax^(t, t^)) and the mean displacement 
(Ax(t,tu;)) denotes an ensemble average over the tracers and over the initial time to which 
is varied over a small time interval centered at t^. We choose a small coarse-grained time 
interval compared to t^, so that we can ignore the aging effect in such short time. We then 
regard all the trajectories in this interval region as having the same waiting time tw In 
practice, at a given t^ and a given At = t — t^ = ti — to, (t > t^,ti > to), we use the 
common method (Q]' pages 118-122) of opening a small time window [t^ — t^, t^ -|- t^] and 
perform an average over all the At with to G [tw—U, +tr\ and ti G [t^ — t^, tmax] in order to 
measure the transport coefficients. We note that t^ < t^ and that t^ax is the maximum time 
of our measurements. The diffusion and mobility constants are then obtained by averaging 
not only over the tracers but also over the initial time to. This common technique allows us 
to obtain an estimation of the diffusion constant and mobility constant by using 82 tracers 
and 75 tracers respectively for this particular system. We check that by reducing the region 
size tr we obtain the same behavior of D{t^) and M{tyj), so that the transport coefficients 
are independent of the averaging technique. 

We collect all the tracers' trajectories from several measurements since our system is well 
reproducible as described in Appendix ^ In each measurement, we consider tracers that 
are separated by at least 10 particle diameters from the sample boundary and from each 
other. This guarantees that the tracer-tracer magnetic interaction can be ignored. 

APPENDIX D: At-LINEAR REGIME 

In this appendix we address the issue of the proper determination of the long At-linear 
regime in the MSD where the diffusivity is calculated. The crossover from the cage dynamics 
regime to the long time diffusive regime seems to be very gradual in Fig. ^ (see for instance 
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FIG. 8: Mean square displacement of tracers as a function of At for the sample A at = 3000s 
and tw = 6000s in (a) log-log plot and (b) linear-linear plot. The dashed lines are linear fittings of 
the At-linear regime. 

the case of = 6000s). We notice that due to the existence of the plateau in the MSD, the 
fitting to the long At-linear regime is of the form (Ax^) ~ A(t^) + 2D{t^)At, where A{t^) 
is a constant dependent only on and D(ty^) is the diffusion constant. Thus, in the log- log 
plot of Fig. d such a linear relation would not be a complete straight line, [note that the 
logarithmic plot is necessary to observe the different time scales in the experiment, and the 
fitting constant v4(t^) is the y-intercept of the fitting of the long At-linear regime, for which 
there is no particular physical meaning]. 

In order to address the issue of the existence of the long-time linear regime we compare 
a log-log plot of the MSD with a linear-linear plot of the same quantity for the sample A in 
Fig. IHK and Fig. |Hb, respectively. We check our results for ty, = 3000s and = 6000s. 

First we note that the log-log plot of Fig. naturally emphasizes the short time scales of 
the cage dynamics. Thus, a linear fitting function plus a constant, (Ax^) = A{t^)+2D{tw)At 
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would be represented as a smoothly varying curve for long At in such a log-log plot, giving 
the impression that the long-time linear regime is not well-defined. In order to determine 
more clearly the large time regime we have plotted the MSD in a linear-linear plot in Fig. ^jp. 
The At-linear regime and the fitting region can be more clearly seen from the linear-linear 
plot. The smaller are also plotted in the figure (t^ = 3000s). As discussed above the 
linear-linear plot shows more clearly the existence of the linear regime while the log-log plot 
shows the existence of the cage dynamics regime. 

We also note that the fast motion inside the cage is determined by the viscosity of the 
surrounding fiuid, and not by the macroscopic viscosity which determines the Tef± at long 
times. Indeed, the motion inside the cage is supposed to be equilibrated at room temperature 
and satisfies the regular FDT according to previous studies, for instance see simulations of 
Lennard- Jones liquids 2l|. The time scales for this microscopic viscosity to dominate is of 



the order of 10~^s, as obtained in light scattering experiments of colloidal glasses |25|. These 
time scales are too fast to be studied with our visualization schemes. Our measurements 
start at At ~ Is and go on up to lO^s. In this long time regime the microscopic viscosity does 
not dominate the dynamics and indeed we find the Tgff is dominated by the macroscopic 
viscosity given by the particles. 

From Stokes-Einstein relation, F = GiramTjv, with F ~ lO^^^A^, rj ~ 10~^N.s.m~^ the 
viscosity of the decalin-bromide solution, ~ 3/im the tracer's diameter and At ~ 10~^s, 
we find Ax ~ 0.002/xm. At this length scale the particle experience the viscosity of the 
surrounding liquid alone. Since we measure much larger displacements, we conclude that the 
tracers are experiencing the macroscopic viscosity and the effective temperature is therefore 
well defined. 

Therefore, we believe that we have achieved the necessary separation of length scales 
which allows us to define a T^tf. The small displacements observed in the response function 
are due to the fact that the magnetic forces ought to be small enough to achieve the linear 
response regime. The small displacements obtained in our measurements are consistent with 
previous confocal microscopy studies of colloidal suspensions near the glass transition j^^, 
which found motions of the order of 1/10 of the particle diameter. 
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FIG. 9: Parametric plot of correlation versus integrated response function for dilute sample C. 
The fitting line shows that the temperature of sample C, T = (323 it 30)K, is close to the bath 
temperature (Tbath = 297K) in our experiment. 

APPENDIX E: STUDY OF SAMPLES C 

Figure El shows the parametric plot of correlation function C versus integrated response 
function ksX for the dilute sample C, (f)c = 0.13. We find M = 3.28 x lO'^ms-^N-^ and D = 
1.46 X 10~^'^m^s~^, corroborating that this sample is equilibrated at the bath temperature 
T = (323 ± 30)K within the uncertainty of the experiment. 



APPENDIX F: LINEAR RESPONSE REGIME 



It is important to test for a well-defined linear response regime for small enough external 
forces where the mobility becomes independent of the force. Indeed, previous work Q did 
not find a linear regime, but instead there is a threshold force below which the particles are 
trapped even under the influence of the external force. On the contrary, here we find that if 
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the observable time is larger than the cage life time, particles are always able to explore the 
structural relaxation leading to cage rearrangements, either in the presence or in the absence 
of the external force and therefore the linear response regime ensues. This is because the 
present experiments focus on larger time scales than previous work. 

In order to find the linear response regime, we apply different forces for samples A and B. 
The results are shown in Fig. ITUk and Fig. ITUb . The collapsing of the results onto a single 
curve of {Ax)/F suggest a linear response for lower forces. The insets show the relationship 
between the velocity, = (Ax) / At, of the tracers along the force direction and the force F, 
which indicates that the external force cannot be larger than 4 x 10~^^N and 6 x 10~^'^N to 
keep the linear response for samples A and B, respectively. It then confirms that the force 
F = 1.7 X 10~^^N used in our studies is in the linear response regime. 

It is important to investigate the age-dependence of the linear response. As we addressed 
in the original manuscript, many measurements have to be taken to reveal the aging behavior 
for one value of the external force. Therefore to measure the aging effect in Fig. we 
measured over 75 trajectories to achieve a statistically reliable average. In order to find the 
region of forces in the linear response regime, many experiments have to be run at different 
forces. Thus, for the determination of the linear response we are obliged to use a smaller 
number of tracers (roughly 10). In this case, in order to improve the statistical average, we 
average over the waiting time and obtain a single velocity for each force, as shown in 
Fig. ^1 (The aging dependence is ignored, and the response functions as shown in Fig. ^1 
are calculated by regular time-translation-invariance methods). This average is only applied 
to find the value of the external force in the linear response regime. Once this force is 
identified, then we perform the full aging analysis by measuring a larger amount of tracers. 

APPENDIX G: SEDIMENTATION OF MAGNETIC TRACERS 

The density of the magnetic tracers is very close to the background PMMA particles. 
Since the sedimentation is proportional to the mismatch between the density of the tracers 
and the surrounding liquid (which is the same as the PMMA density), then we expect that 
the effect of sedimentation would be small. 

In principle, this effect is negligible in the more dense sample B since the tracers remain 
in the x-y plane for few hours without any noticeable sedimentation effect. In the case 
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FIG. 10: The response functions for (a) sample A and (b) sample B. The insets are plots of the 
velocity of magnetic beads as a function of the force, Vx versus F. The straight line indicates the 
linear response regime. The velocity error is due to the accuracy to locate the bead positions, 
smaller than the symbol size in the inset of (a). 

of the less concentrated sample A, sedimentation of tracers can be seen after one hour of 
observation. Since we measure up to = 6000s, it is important to check the effect of 
sedimentation on the diffusion of particles for this sample. 

Unfortunately, we cannot test directly the effect of this sedimentation, since we cannot 
match the density of the tracers with that of the fluid and PMMA particles (the tracers 
are magnetic and therefore slightly heavier than the PMMA particles by construction). 
However we can indirectly test the main issue at hand: What is the effect of an external 
driving constant force of the order of gravity (~pN) on the tracer diffusion? 

To this end we apply a magnetic external force of the same value as the sedimentation 
force and measure the resulting effect on the diffusion constant for sample A. We found that 
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At (s) 



FIG. 11: MSD for the sample A at (a) = 500s and (b) = 6000s analyzed from the data 
without magnetic force and with magnetic force, respectively. 

the effect of the external force on diffusion is negligible as long as the force is small enough 
(~pN). Figure HU plots the MSD for the sample A at (a) = 500s and (b) = 6000s. 
We obtain the MSD analyzing the trajectories of tracers without magnetic force and with 
magnetic force (in this case, we properly subtract the average value of the displacements). 
The figure indicates that the diffusion is the same with or without the external force, thus 
confirming its negligible effect for this case. 

APPENDIX H: DETERMINATION OF THE AGING EXPONENTS 

As described in the main manuscript, all the MSD for different aging times collapse 
onto a master curve when plotted as f^Cit^.t^ + At) versus AtjtJ^ . In order to search 
for the best value of the aging exponents a and (5 that give the correct data collapse, we 
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FIG. 12: Log plot of standard deviation, In cr^, which quantifies the cohapsing of t'^C{t^,tu, + At) 
for all the data of sample A (in color scale or gray scale) . The minimum of cj^ (black basin) gives 
the best aging exponents at a + /? = 0.34 ± 0.05, (3 = 0.48 ± 0.05. 

calculate the standard deviation of all the curves as defined below: 

E E 

ln{At/ti)Htn,) 

[ln(CC(t., At/t^)) - ^ E HtZC{t^,At/t^J)r. 

Htw) 

We search a wide range of a and f3 to minimize a^. The result is plotted in Fig. IT^ 
showing that the best data collapsing with = 0.1 can be achieved ai a + j3 = 0.34 ± 0.05 
and f3 = 0.48 ± 0.05 as reported in the main text. 

Spin-glass models predict a general scaling form Cag(t,t«,) = Cag{h(t)/h(tu_,)), where h{t) 
is a generic monotonic function Qj. In particular, a familiar form is h{t) ~ where fi is 
a power law exponent giving rise to the cases of simple aging, /i = 1, super aging, /i > 1, 
and sub-aging, /i < 1. Assuming the above form, the linear behavior Cs^gitw + At,t^) ~ At 
for larger At requires fi = 1. This is just a special case of our scaling ansatz Eq. ()5aj) with 
P = 1. As we show in Fig. [T21 our data cannot be collapsed for this parameter. 
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APPENDIX I: SCALING BEHAVIOR OF LOCAL FLUCTUATIONS 



The local correlation function C and local integrated response function x studied in the 
manuscript are calculated from each individual tracer trajectory. In a condensed colloidal 
sample, the tracer trajectories are always confined at a local position. Therefore, the cor- 
relation C and response x fo^' individual particle can be regarded as the coarse-grained 
local fluctuations of the observables as investigated in 31 1. 

Furthermore, in order to improve the statistics of our results, the PDF P{C) is calculated 
not only for all the tracers, but also over a time interval much smaller than the age of the 
system. In practice, we calculate the i—th tracer's local correlation function as Ci{t^ + 
At,t^). At a given At, we open a small time window [At — ts,At + tg] [tg ^ At) and 
count all the Ci(t^ + At-,-, t^) with Atj G [At — t^. At + ts] into the statistics of P{C). The 
calculated P{C) is a mixture of the local and the temporal fluctuations of the observables. 
Similar technique is performed to calculate P{x)- 

Below we derive the scaling law for the PDF of the local correlation function. Let us first 
recall the scaling behavior of the global correlation function in Eq. @, 



tlC{t^ + At, t^) = /,,(At/t^^), (II) 

where we add a bar to C to distinguish the global correlations from the local C. The average 
is taken over all the tracer particles. We can rewrite the global correlation function as the 
integration of the PDF P{C) as: 



C = j P{C)CdC. (12) 
Furthermore, we obtain the relation of //) and P{C) by substituting Eq. (|I2|) into Eq. 



fD{At/tJ) =tZC = J t-^P{C){tZC)d{tZC). (13) 

For a given At/t^^, is equal to a constant, and Eq. (jT3)) requires that t^°'P{C) should 
only depend on t°C. In other words, t^'^PiC) is a function of At/tw^ and t^C. Then we 
define Fd as: 

Fn{tlCAt/ti)=t-JP{C). (14) 
31 



A similar formula can be obtained for P{x)- Eventually, we obtain the scaling ansatz of 
P{C) and P{x) shown in Fig. |3J 

P{C)=tZFn{tZC,At/ti), (I5a) 

Pix)=tZFMKx,^t/ti), (I5b) 

where the universal functions FD{x,y) and FM{x,y) satisfy 

j FD{x,y)dx = foiy), (I6a) 

Fm{x, y)dx = fuiy). (I6b) 



APPENDIX J: STUDY OF THE PDF OF THE AUTOCORRELATION FUNC- 
TION 

The PDF of the autocorrelation function follows a modified power law t^"P(C) oc (t^C + 
Cq)~^, as we see in Figs. H^, and|3t. From Eq. (jTSl), we have 

tic = j^^^\~'^P{C){tlC)d{tlC) 

_ j^'-\x + CoY'-xdx ^"^^^ 
" j^^-\x + Co)~^dx' 

The cutoff Ccut {Ccut ^ C*o) is introduced to make the integral converge and we always take 



A > 1, then 



CC^=^[l-(A-l)(^ + ir^]. (J2) 
For A > 2, the last term in Eq. (p2|) is negligible and t^C ~ Co /(A — 2) mainly depends on 



the short-time parameter Cq. For A < 2, we have t^C ~ |ztCo(%^)^ ^ and the long-time 



parameter Ccut dominates. 

Following the previous discussion of P{C), we define the i-th tracer's local correlation 
function as: 

Ci{t^ + At, t^) = {Axit^ + At, t^)'),/2, (J3) 

where the average {■■■){ is calculated for only the i-th tracer's trajectory by opening a small 
time window [At — t^. At + 1^] [tg ^ At) at a given At, and counting all the Cj(t^ + Atj,t^) 
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with Atj G [At — ts, At + ts] into the statistics of P{C). We should note that if there is no 
average P{C) can be reduced to a simple form: 

P{C) = P((Ax2),/2) P(AxV2). (J4) 

This form can be further reduced since 

P(AxV2)rf(AxV2) = P{Ax)dx, (J5) 

therefore at the limit of the average 0, 

P{C) = P{{Ax^)i/2) P(AxV2) = P(Aa;)/Ax. (J6) 

Therefore our P{C) is relate to the probability P(Ax) usually studied in previous works 



24| . We find that P{Ax) can be approximately fitted by a broad tail power law behavior 



for large Ax, consistent with equation [J6l However the data of P(Ax) can be better fitted 



by stretch exponential rather than power law, which is consistent with previous works |24 1 . 
On the contrary, P((Ax^)i/2) can not be fitted by stretch exponential. Thus we believe 
that the power law fit is more proper to describe our result of P(C). The difference with 
other studies might be due to the fact that our samples are glassy while others work in the 
supercooled regime. 
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